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Abstract 

In QCD hard scattering cross sections, the color content of the underlying 
hard scattering evolves with a factorization scale. This evolution is controlled 
by an anomalous dimension matrix, specific to each hard-scattering reaction. 
Anomalous dimensions are determined from the renormalization of products of 
ordered exponentials of the gauge field, which describe the coherent radiation of 
gluons by incoming hadrons and the observed jets or particles of the final state. 
The anomalous dimensions depend on the kinematics of the underlying hard 
scattering, but are free of collinear singularities. A number of these matrices 
are available in the literature. Here, we exhibit one-loop mixing matrices for 
the full list of 2 — » 2 reactions involving light quarks and gluons. The eight-by- 
eight anomalous dimension matrix for gluon-gluon scattering shows a simplified 
structure in the basis corresponding to definite color exchange in the t-channel. 
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1 Introduction 



Inclusive, short-distance hadronic cross sections may be factorized into convolutions 
of universal, nonperturbative parton distributions and/or fragmentation functions, 
and perturbative hard-scattering coefficient functions. In this procedure, a factor- 
ization scale fi separates nominally long- and short- distance components of the cross 
section. When /i is chosen large enough that the coupling a s {jj?) is perturbative, 
it is possible to calculate the coefficient functions, and hence to determine the cross 
section, assuming that the relevant nonperturbative densities are known. 

Factorization is an expression of quantum mechanical incoherence. The parton 
distributions summarize the internal dynamics of the incoming hadrons, and of their 
fragments that evolve into the final state. These distributions are independent of 
each other and of the fragmentation of scattered partons into jets and hadrons, in the 
final state. They are also incoherent with the short- distance dynamics that produces 
the hard scattering. Finally, they are incoherent with the dynamics of "noncollinear" 
gluons of any energy - from the hard scale Q down to zero - that are radiated away 
from the incoming hadrons, or from the observed components of the final state Jl], ^ 

Each of these components of a hard-scattering process: parton distributions, par- 
ton (or jet) fragmentation, hard scattering and soft gluons, may be thought of as 
governed by an effective field theory, in which the parton dynamics at differing scales 
and directions has been removed, and sometimes replaced by specific operators. At 
fixed order, coefficient functions combine contributions of quanta that are off-shell by 
the order of Q with those from noncollinear soft gluons, because divergent soft-gluon 
behavior cancels, leaving finite contributions at each order in perturbation theory. 
These contributions, however, retain sensitivity to momentum scales much smaller 
than fi, and may produce large, although finite, higher-order corrections. 

Noncollinear soft-gluon dynamics is an important ingredient in a variety of re- 
summed cross sections. Their first complete treatment was given by Collins and 
Soper for the transverse momentum (Qt) distribution of hadron pairs in e + e~ anni- 
hilation 0]. Because this cross section is not fully inclusive, it is necessary to resum 
(Sudakov) logarithms of Qt due to soft gluon emission to determine the cross section 
near Q T = 0, a procedure intensively studied for Drell-Yan pairs [§, |J. As is typical of 
such resummations, nonperturbative effects remain important phenomenologically at 
realistic energy scales. In fact, the perturbative resummation suggests the functional 
form of these effects 0. A similar resummation is necessary to define elastic scatter- 
ing amplitudes for hadrons and partons at large momentum transfer and fixed angle 
|0, |, H, |T0|| . In both of these cases, leading logarithms are generated by soft gluons 
emitted collinear to incoming partons. Nevertheless, a treatment of noncollinear soft 
emission, i.e. soft gluons emitted at angles between the directions of the primary 
incoming and outgoing hard partons, is necessary to extend the formalism to next- 
to-leading logarithm and beyond. The pattern of this emission is essentially coherent 
TT|; as we shall see below, it is determined primarily by interference between emis- 



sion from different hard partons. Despite this, it remains incoherent with the collinear 
dynamics of incoming and outgoing hard partons. That is, it factorizes from parton 
distributions and from suitably defined final-state jets. 
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In inclusive quantities such as the Drell-Yan cross section da/dQ 2 , soft gluon 
dynamics is somewhat hidden in higher order corrections to short-distance coefficient 
functions. It appears indirectly, through integrable logarithms at the "edge of phase 
space" , where the observed final state uses up all partonic energy, leaving none over 
for QCD radiation. These phenomena have been studied in inclusive electroweak 



reactions such as the Drell-Yan cross section |12| , |13| , and more recently in inclusive 
QCD reactions such as heavy quark [14, 15, 16, 17, IS, 19, 20] and dijet production 

Beyond leading logarithms, all of these resummations require a treatment of soft, 
noncollinear, gluons. The hard scattering of partons of specific flavors involves, in 
general, a number of possible color exchanges. On quite general grounds, we expect 



this information to be reflected in observable properties of the final state p2[. At the 



same time, it is clear that if soft gluons are factorized from the hard scattering, the 
color content of the latter depends on the scale at which the separation is carried out. 
Nevertheless, physical cross sections cannot depend on this scale. The combination of 
factorization, and invariance under scale choice ensures the possibility of resummation 



in the scale dependence [23 



In this procedure, the emission of noncollinear soft gluons is described by matrix 
elements of products of ordered exponentials of the gluon field, coupled at a point 
representing the hard scattering [|7|, [17], [l^, ^(J, |2l| . The factorization scale between 
the hard and soft gluon functions may be interpreted as the renormalization scale for 
these operators. If the physical cross section is defined so that the soft gluon function 
depends on a single ratio of a phenomenological mass scale to the renormalization 
scale, then the entire soft function is determined by its anomalous dimensions. 

Rather than review the specifics of QCD resummations in detail, we shall treat 
here one feature common to them all, the evolution of color exchange in noncollinear 
soft gluon emission. Our discussion will therefore treat gluons interacting with ordered 
exponentials instead of physical partons. We will present the full set of anomalous 
dimensions relevant to the scattering of light partons, including quark-quark, quark- 
antiquark (which have been presented previously), quark-gluon, and (by far the most 
complex case) gluon-gluon. Each of these partons will be replaced, for this discussion, 
by Wilson lines of the corresponding color representation. The justification for this 
approximation, and the detailed role that it plays in each cross section, may be found 
in the references cited above. 

We begin, in Sec. 2, with a discussion of the renormalization of products of ordered 
exponentials, with the aim of isolating the evolution of color mixing. We shall see 
that, although the analysis immediately encounters collinear divergences, it is easy 
to identify finite anomalous dimensions that govern the mixing of color. We go 
on, in Sec. 3, to recall, for definiteness, the soft gluon contribution to inclusive jet 



cross sections in hadronic scattering pl| . This "soft function" is manifestly free of 
collinear singularities, and its evolution, with mixing, reflects the color content of 
the underlying hard scattering. In Sec. 4, we introduce a color basis notation that 
is convenient for the solution of the evolution equation for the soft function. This is 
followed, in Sees. 5 and 6, by the computation of the anomalous dimension matrices 
for various processes, along with their eigenvalues and eigenvectors. We have included 
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two appendices that describe some of the details of the calculations. 



2 Soft Anomalous Dimensions from Wilson Lines 

Our discussion begins with ordered exponentials, or Wilson lines. A general Wilson 
line is labelled by a path in space-time, C, beginning at point z and ending at point 

z' 



W[C;z',z] = Pexp 



-ig f X2 d V ^--A(y(rj)) 



X 



The variable of integration r] parameterizes C, with y(X\) = z, y(\2) = z' . 

Wilson lines have myriad applications in quantum field theory, but their interest 
here is as a model for the radiation of soft gluons by fast-moving partons (quarks 
and gluons). Any hard-scattering event, be it based on electroweak or QCD hard 
scattering, results in nonabelian radiation, in much the same manner as in classical 
electrodynamics. Indeed, as in classical electrodynamics, the pattern of soft radiation 
is determined by the charge currents at long times before and after the scattering 
event. The factorization properties of inclusive hard-scattering cross sections in QCD 
rely on this feature 0. 

Initial and final-state partons are represented by ordered exponentials in which C 
is a straight line, in the direction of the four-velocity, /3, of the relevant parton. For 
inclusive cross sections, the lines extend from a common point x to infinity, either 
from the distant past, or toward the distant future. We shall adopt the notation of 



Ref. [^T] for such Wilson lines, 

$ { /\\ 2 ,\ 1 ;x) = Pexp(-ig d V (3-A^( v (3 + x)\ , (2) 



where the gauge field A^> is a matrix in the representation of the flavor /. Comparing 
this definition with Eq. (fj), we have y{rj) = rj/3 + x. 

Although the Wilson lines $ in Eq. (0) extend to infinity, they reproduce the dy- 
namics of the gauge field for only a finite time before and after a hard scattering. At 
long times, nonperturbative dynamics dominate, nonabelian charges are neutralized, 
and a picture in terms of isolated currents fails utterly. Therefore, the procedure ap- 
plies only to cross sections for which such long-distance effects cancel in perturbation 
theory. The final answers for such quantities are insensitive to the infrared properties 
of amplitudes generated from Wilson lines of the form of Eq. (0), although it may be 
necessary to deal with these properties at intermediate steps in any calculation. In 
fact, it is their ultraviolet behavior that we shall find most interesting. 

This approach has been applied to the classic electroweak hard scattering pro- 
cesses, e + e~ annihilation, DIS and Drell-Yan vector boson production |L2|, [13], j24"f . 
In each of these cases, the relevant operator is a product of two $'s coupled at a 
color singlet vertex. Equivalently it is a single Wilson line, whose direction changes 
discontinuously at a single point. For e + e~ annihilation, the curve C consists of two 
straight paths, both extending to positive infinity in time, describing the creation of 
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a quark pair at a point. For DIS, C comes in from negative infinity and changes 
direction, describing the scattering of a quark or antiquark, while for Drell-Yan, C 
represents the annihilation of a pair. 

For definiteness, we consider the Wilson line analog of the Drell-Yan cross section 



24[, where the configuration just described corresponds to the product 

<?W = S aia2 ^{0,-oo;x) a2b2 ^{0,-oo;x) aibl . (3) 

The vertex at x is associated with ultraviolet divergences, as are self-energies of the 
eikonal lines. The divergent one-loop diagrams are shown in Fig. [I]. These vertex 
and self-energy diagrams are conveniently renormalized by perturbative counterterms, 



corresponding to multiplicative renormalization of the composite operators |25|, ^6 



W™°\x) = ( a s , ^Jl) W^\x) , (4) 



where jy( DY )(°) is the bare operator, and where we have used the invariance of W 
under rescalings of the velocities In the following, we suppress dependence on 
x, the position of the composite vertex. The independence of the unrenormalized 
operator from the renormalization scale /i implies that the renormalized operator 
satisfies the renormalization group (RG) equation^, 

^^ = rJ« s M] w m). (5) 
dp v MS 



The anomalous dimension Yw{ot s , A) is the same for all color singlet products of the 
form of Eq. ([|), and is sometimes referred to as r cusp p7| . 

Because W^ DY \a s , Pi, n) is a gauge-invariant operator, so is Y w . A direct calcula- 
tion verifies that, as suggested by Eq. (|5|), Y w is somewhat unusual for an anomalous 
dimension, in that it diverges in the limit of light-like velocities for the external lines, 
j3f = 0. Fixing the f3f to be small, but nonvanishing, we find the explicit expression 



W — 

7T 



In 



-2/3! • p 2 



(6) 



The divergence at light-like velocities is a collinear enhancement, due to the emis- 
sion of gluons in the direction of the eikonal line. In suitably-defined inclusive cross 
sections, the corresponding collinear divergences either cancel, or are absorbed into 
parton distributions 0. We shall see that the same is true for a product of two Wil- 
son lines like Eq. (§) and its generalizations. Specifically, the collinear divergences in 

1 Strictly speaking, the multiplicative renormalization in Eq. (^), and the corresponding RG equa- 
tion (^|) apply to the renormalization of the composite operator only; the renormalization associated 
with the QCD Lagrangian is assumed to have been carried out. In the renormalization proof 
was formulated in terms of the vacuum expectation value of the corresponding loop. 
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IV are universal, and may be factored from the underlying local coupling. In com- 
puting loop corrections to W and $, we treat the eikonal lines as on-shell particles, 
normalized by the square root of the residue of the eikonal two-point function. This 
ensures gauge invar iance. 

To demonstrate factorizability, we shall find it useful to work in an axial gauge, 
n ■ A = 0, with n • fa ^ 0, i — 1, 2. In any gauge, a single Wilson line $i in quark 
representation, defined on a ray in the direction of velocity (3 is also multiplicatively 



renormalizable 25 



(e)(0) 




(7) 



where again, rescaling invariance in the velocity and the gauge vector limits possible 
arguments of Z. Eq. (0) implies that $^ satisfies a renormalization group equation 



d 



(?) 




(?) 



with another anomalous dimension. Clearly, requires the same renormaliza- 



tion constant. 

The one-loop expression for T$ in axial gauges has the same collinear (/3 2 
divergence as IV, but now through gauge-dependent logarithms, 



a. 



1 1 -n[ 

(3 2 \n 2 



— Of 

7T 



In 



2\(3-n\ 
/(3 2 \n 2 



0) 



(9) 



Evidently, the collinear divergences in this anomalous dimension exactly match those 
of r^- 1 , with the opposite sign. As we shall see, this result is quite general. 

To make use of this correspondence, let us define a new operator, which generates 
the sum of diagrams, containing the Drell-Yan vertex, that cannot be reduced by cut- 
ting a single eikonal line. We do this by simply dividing out the vacuum expectation 
values of the two external Wilson lines in Eq. (|3|), which defines W^ DY \ 

S (DY) = W (BY) ( |$jf>(0, OIO)- 1 (0|*g(0, -oo; OJIO)- 1 , (10) 

where in the vacuum expectation value a trace and an average of color of the operator 
$ is assumed. 

Given the RG equations for W^ DY ^ and the $'s, the function 5 ,( - DY ^ satisfies a 
renormalization group equation of exactly the same sort as W( BY \ 

d 



dfx 



S^=T s (a s ,fa,n) 



(11) 



with a "soft" anomalous dimension, 



w 



r (?) 



7T 



O'.s 



In 



r (9) 



O's 



(3 2 \n 

-fa -fa \n 2 \ 
2\fa ■ n\ \fa ■ n\ 




(12) 
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We may think of T$ as the "coherent" part of the "cusp" anomalous dimension T^, 
that is, the part that depends on the specific color-singlet vertex in Eq. @, but 
is independent of the mutually incoherent collinear behavior of the lines themselves. 
Our definition of T$ is, however, somewhat arbitrary, and we might change it through 
a finite renormalization and/or a change of gauge. A natural choice is A = gauge 
in a frame where /3f = 5 M+ and (3% = <V-' m which case Ts = (a s /ii)CF (1 — «vr) fll2| 



but we are even free to make T^ Y ^ vanish. Any choice leads to the same amplitude. 

The value of this approach is found by generalizing W/( DY ) to products of multiple 
Wilson lines. For QCD hard-scattering processes, products of four lines are most 
interesting, and illustrate the general case. As noted above, the collinear divergences 
represented by Eq. (^) are more general than the color-singlet vertex in W^ DY \ pre- 
cisely because they may be factored from this vertex, and from each other. In the 
anomalous dimensions, which appear eventually in exponents, this factorization is 
guaranteed to be additive. 



In the context of QCD hard-scattering |17| [18], [K], g0| |21| , most of the above 
considerations apply, although now with a more complex and interesting color flow. 
In this paper, we shall be concerned specifically with soft radiation for cross sections 
due to 2 — * 2 scattering of light quarks and gluons at short distances, of relevance 
to high-pr jet production. Analogous reactions are also at the basis of heavy quark 
production. In this larger context, the value of isolating single-logarithmic soft gluon 
dynamics from corresponding, universal, collinear dynamics is clear. We shall review 
very briefly in Sec. 3 the role of products of Wilson lines in hard cross sections, using 
as an example the resummation of threshold corrections in jet cross sections. First, 
however, we show how the generalization of the color-singlet product W^ DY ^ leads 
naturally to a matrix of anomalous dimensions that describe the mixing of color with 
changing scales. 



Following Ref. [pijl , we introduce a class of operators coupling four Wilson lines, 



corresponding to the primary partons involved in a hard scattering, 

(13) 

and a corresponding irreducible soft function, 

S? = WP J] (01^(0, -oc^lO)- 1 nWSfWjOJlO)- 1 . (14) 

i=A,B i=l,2 

As above, represents the four-velocities of initial- and final-state partons. The 
superscript f represents a scattering process, 

f A {£ A , r A ) + f B (i B , r B ) -> /i(pi, n) + f 2 (p2, r 2 ) , (15) 

involving four partons, with flavors L, incoming momenta Zj and outgoing momenta 
Pi, and color n. The color tensor (cf^ describes the couplings of the ordered 
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exponentials with each other in color space. If, for example, the incoming lines are 
quark and antiquark, the indices r iy dj are in the fundamental representations of 
SU(3). Examples of bases for the color tensors cj then include: singlet and octet 
matrices in the s- or in the t-channel, or s-channel color singlet and t-channel color 
singlet. 

The relationship between the Wj and the Si is shown in Fig. |2| in axial gauges. 
The left-hand side of the figure shows a general momentum configuration of a "leading 
region" that contributes logarithmic behavior to (the vacuum expectation value 
of) Wp; lines have been grouped according to their momentum flow. Such configura- 
tions may be identified by standard Minkowski-space power counting |§ . As indicated 
in the diagram, all lines fall either into "jet-like" subdiagrams, labelled 1, 2, A, B, or 
into a "soft" subdiagram, labelled S. The jet subdiagrams consist of lines that are 
all on-shell with light-like momenta parallel to the directions of the corresponding 
The momenta of lines in S vanish in all four components. As indicated in the figure, 
in each such "leading" region, all collinear contributions factor into the self energies, 
in axial gauges. Demonstrations of this factorization, based on Ward identities in 
QCD may be found in Refs. H, M and M in the case of jets of partons. The eikonal 
lines in Fig. |2] obey the same Ward identities as do quark lines, and do not change 
these arguments. 

({) a) 
Under renormalization, the set of operators Wj , and hence the S} , will mix in 

general. Because the eikonal lines do not include recoil effects, they cannot annihilate, 
and there is no mixing between vertices that link other sets of Wilson lines (corre- 
sponding to other flavor flows f). The corresponding renormalization group equation 
is now in terms of a matrix of anomalous dimensions, 

v-j-wP = (rU) ;; KA,ft) wf ] 

„ d o(f) (r (t)\ ( n Pi-fa \ n2 \ \ o(f) rifi N 

dfi V \Pi ' n \\Pj ' n \J 

where in the arguments of Fg, i ^ j. The two anomalous dimension matrices are 
related by a generalization of Eq. ( jl2[) for Drell-Yan, 

rV(a s , Pi, n) = Y%j {a s , (3 h ft) - S u ^ Tf } I a s , ^L^L ) . (17) 



The anomalous dimensions for Wj include collinear divergences in the limit of van- 
ishing (3f. These divergences, however, are universal, and may be absorbed into the 
renormalization of the $'s. This is possible precisely because any collinear divergence 
may be absorbed into a self-interaction of one of the Wilson lines, as argued above. 
Since the renormalization of self-energy graphs does not mix colors at the vertex, 
collinear divergences in , Eq. ([TED , can only appear times the identity in the space 
of color tensors. Therefore the soft anomalous dimension matrix Tg jj summarizes 
the mixing of operators under renormalization, but is free from collinear divergences, 
as long as /3j • f3j ^ for i ^ j. The soft anomalous dimension will be, as in the 
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case of r^ Y ^ above, free of universal collinear divergences, not only at one loop, but 
to all orders in perturbation theory. 

Note that the choice of finite renormalization for the $'s, and of overall gauge, 
encountered above, may be interpreted as a renormalization scheme for defining the 
composite vertex and its soft anomalous dimension. In the next section, we shall 
exhibit a class of infrared safe soft functions that occur in the threshold resummation 
of dijet cross sections [2I|, and whose evolution is controlled by the tensors Tgjj- 



3 Soft Functions in Dijet Production 

In this section we recall the role that eikonal operators like those discussed above play 
in the resummation of threshold corrections in dijet production, 

h A {p A ) + h B {p B )^Ji{Pi) + J 2 {p 2 ) + X(k) , (18) 

at fixed dijet invariant mass, 

Mjj = (pi + p 2 ) 2 , (19) 

and at fixed rapidity difference, 

A, = i In = (») 



with t and u the usual Mandelstam variables, defined as 

* = (PA-Plf 



u = (PA-P2) 2 . (21) 



As discussed in Ref. ETJ, many other choices of cross section involve essentially the 



same physics. Again to be specific, we assume that the two final-state jets in (0) are 
identified by cones of opening angles S^. 

We shall be concerned with the contribution to this cross section from an underly- 
ing partonic scattering like Eq. (|15D . In particular, we shall review some relevant fea- 
tures of the resummation, in this cross section, of "threshold" corrections to the hard- 
scattering function. These are corrections of the general form [ln n (l — z)/(l — z)] + , 
with z = Mjj/s, where s is the invariant mass squared of the partons that initiate 
the process. Referring to Eq. (|1^), we take £ A = %aPa, £b = xbPb, with p A and 
Pb the momenta of the incoming hadrons and x a ,b momentum fractions, so that 
s = 2x a xbP A - Pb- We represent the coupling of soft gluons to the hard partons of the 
flavors f A , fs, fi and / 2 by Wilson lines, as in Eq. (H), tied together at the vertices 
of Eq. (0). In terms of the operators we define a dimensionless "eikonal cross 
section" , describing the emission of gluons by hard partons in eikonal approximation, 



(f,eik) / wMjj 
a LI 

V ^ 



Ay,a s (/i 2 ),^e) = £ 5 ( w -^ Cik) (£,^) 



x(o|f [(wf } (o)) \ bi} ] \0W [wP(o) {bt} ] |o> , 

(22) 
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where £ designates a set of intermediate states. The eikonal cross section depends 
on the rapidity difference Ay (or equivalently, the partonic center-of-mass scattering 
angle) through the relative directions of the Wilson lines coupled at each operator 
We choose the dimensionless weight w in Eq. (]22f) to reproduce the phase space 
available to soft gluons near partonic threshold (z = 1) ||21|| . For the cross section in 
question, this requires that we fix the total energy emitted by the eikonal lines into 
the final state, but outside the cones Si. Denoting this energy as k^, we have 

2k° 

w(eik)(e '^ ) = M^- (23) 

To handle the collinear divergences of cr^ f ,clk -* due to initial state interactions, and to 
avoid infinite energy emitted into the cones, we may think of the eikonals as having 
small "mass", /3f. We shall not need this regulator in any of our explicit calculations 
below, however. 

As in the previous section, we want to isolate the contributions of noncollinear 
gluons, which control color mixing. For this purpose we factorize the eikonal cross 
section, Eq. (^), in the same manner as a hadronic cross section @, 21|. Initial-state 



collinear divergences can be factorized into eikonal "parton distributions", which we 
label Jin- Analogously, final-state collinear enhancements may be factorized into 
outgoing "jet" functions, j'out- 

The result of this procedure is to write the eikonal cross section as a convolution 

(f) 

of an infrared safe soft function S L j with the functions j'in and j'out > 

(f,eik) / wMjj 2 \ 

I — - — , Ay, a a {(i ),Oi,ej = 

1 

dw^dw^dw^dw^dw'g 5 (w — w[ — w' 2 — w' A — w' B — w' s ) 

c=A,B \ f 1 ) d=l,2 \ A* / 

xS!}(^,AyM» 2 )) ■ (24) 

Here w[ = 2k'j°/Mjj, with kf the energy of partons emitted outside the cones, and 
associated with the function labelled by % — 1, 2, A, B, S. We are using the incoherence 
of these "effective theories" , as mentioned in the introduction. 



For completeness, we give definitions for the incoming eikonal jets |T7L ETJ: 



x (0| Tr{ T[$#>t( , -oo; y)]T[$^(0, -oo; 0)] } |0> , 

(25) 

with i — A,B and y v = (yo, 0) a vector at the spatial origin. The argument e = 4 — D, 

( f ) 

in D dimensions, denotes the presence of collinear singularities in . Similarly, for 
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the outgoing eikonal jets we have ||21|| , 



,(/i)t/~, n. nMlt\/tlTT<T.(/i)i 



x(0|Tr T[$^ (oO;0 ;0)]|0^|T[$^(oo,0;0)] |0> , 



(26) 



with i — 1,2. In this case, collinear singularities cancel due to the jet definition in 
the sum over final states, and the infrared regulator e is replaced by the cone size Si. 

Mellin transforms of the eikonal cross section with respect to the weight conve- 
niently isolate contributions from noncollinear soft gluons that are singular at partonic 
If we take the Mellin transform of Eq. (|3|) , 



threshold [21 



dw (1 



w 



N-l (f,eik) 



wM 



a 



JJ 



LI 



Ay,a s (fj, 2 ),5i 



(27) 



the convolution in the eikonal weights becomes a product of moments, up to cor- 
rections that vanish as a power of the moment variable N for large N. Hence the 
moments of the soft function can be extracted from the eikonal cross section, as 
follows, 

'(f^Ay,^ 2 ),^ 




1 



ln m (l-2) 



1- 



may be reconstructed from logarithms 



All the leading power singularities, 

of the leading power (O(N )) N dependence. 

We recognize Eq. (p8|) generalization of Eq. 

cross section 



[L4D to the moments of the 

By removing the mutually-incoherent jet factors, which in axial gauge 

(f) 

correspond simply to cut parton self-energies, we isolate in the function S^j the inter- 
ference of amplitudes in which gluons are emitted by one eikonal line and absorbed 
by another. In axial gauge, cut diagrams of this topology incorporate what we might 
describe as "interjet" coherent [[111 



radiation in the cross section. 
We have not yet discussed the renormalization of the soft function and eikonal 
jets. In the next section we will see that renormalization leads to the exponentiation 
of the ^-dependence of the soft function, and hence of the energy dependence of its 
transform. 

Now let us relate the soft function Sli to threshold resummation for dijet cross 
sections. To derive singular behavior at partonic threshold, we consider moments 
with respect to r = Mjj/S for perturbative, partonic cross sections, in which the 
external hadrons are replaced by partons in an infrared-regulated theory. In Ref. 
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PH , we showed that, up to next-to-leading logarithms, this cross section factorizes 
under such a Mellin transform, 

(29) 

The factorized cross section is illustrated in cut diagram notation in Fig. [| Here, 
Sli is exactly the matrix of soft eikonal functions constructed above. Hjl absorbs 
quanta, represented by h* L and hj (see also Eq. (|3l|) below), which are off-shell by 
the order of Mjj. The remaining functions if) and J include, respectively, virtual and 
final-state partons parallel to the incoming hadrons (partons) and to the observed 



outgoing jets. Explicit definitions of the jet functions are given in Ref. \TL\. Their 
detailed features are not of interest here, primarily because they are color diagonal, 
and incoherent with the soft radiation included in Sli- Our interest is in the interplay 
of color structures between the soft tensor Sli and the hard color tensor Hjl- 



4 Color Evolution in Jet Cross Sections 

In this section, we exhibit the evolution of color exchange in hard scattering. Although 
our notation is taken from the dijet cross section of the previous section, it is quite 
general. 

4.1 Evolution for Soft and Hard Functions 

The composite operators of the soft function S LI , as we saw before, require renormal- 
ization. The matrices H and S occur in a product in Eq. ( |2"9"| ) and (consistent with 
the discussion of Sec. 2), must therefore renormalize multipicatively, with separate 
renormalization factors for the amplitude and the complex conjugate 0, |23fl , 



=A,B,1,2 



S">S = (Z?<)L B SM; A1 . (30) 

Here Zj is the wavefunction renormalization constant of the iih incoming parton field 
(of flavor /a ■ ■ ■ ^2) in Eq. (]15D, involved in the hard scattering, and Z S CD is a matrix 
of renormalization constants, describing the renormalization of the soft function. 
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In the following, we shall adopt the convention of using Roman indexes to refer to 
the color structures in an arbitrary basis. In jet cross sections near partonic threshold, 
the hard-scattering function is a product of two purely virtual factors [[HJ 



Hfl ^,Ay,a.^)J [^,Ay,a s ^)j hf A,, a s ^) j , (31) 

one from the overall amplitude and one from the complex conjugate. This simplifi- 
cation is due to the fact that near threshold there is insufficient phase space for 
the emission of hard partons, aside from those observed in the hard scattering (and 
hence summarized in the jet functions in Eq. (|29| )). 

The hi in Eq. fl3~lj ) are the coefficients of color tensors in the expansion of the hard 
scattering amplitude hj d2 treated as a matrix in color space, 

/»«.,, A,, = £Ag> p^,A v ,a,(^)) (c») dA A , (32) 

(f) 

with the c K the same color tensors introduced to define the soft functions, through 



Eqs. (|13"D, fl22|) and (^8|). Here we may point out a strong similarity to the effective 
theory formalism, in which the c/s play the role of operators in a "low-energy" theory, 
and the /i/'s are coefficient functions that summarize the high-energy components that 
have been "integrated out" . 

From Eq. fl3"0|), the soft function S^J satisfies the renormalization group equation 

10 

(/% + ^tg) SS = -(^^S - S&Ff)* , (33) 

where we encounter the same soft anomalous dimension matrix, Fp, as in Eq. (|i~6|) , 
which is computed directly from the UV divergences of the soft function. We will 
calculate this matrix for the basic partonic processes f in a minimal subtraction 
renormalization scheme, taking e = euv — 4 — D, with D the number of space-time 
dimensions. The one-loop anomalous dimensions, 

4 f) (^) = -f^Res_o4 f) (^e), (34) 

are obtained from the residues of the UV poles contained in the matrices of renor- 
malization constants of the soft function. 

The moment dependence of the product of the hard and soft functions may now 
be determined directly from Eq. (|33|) [ f2TH , 



T T (f) (m ' A ^ (/i2) J ^(^,Ay,a s (^) 




Tr tf( f ) — ^,Ay,a s (fx 



xPexp 
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X 



SV(l,Ay,a s (M 2 jj/N' 



xPexp 



(35) 



where the symbols P and P refer to path-ordering in the same and in the oppo- 
site sense as the integration variable //. P orders Tg\a s (/j 2 )) to the far right and 
T < s\a s (M 2 J j/N 2 )) to the far left of the products. In Eq. (|35|) the trace is taken in the 
basis of color structures, as 



Tr \H®SV \=H { ?}S { P 



l IL u LI 



i,(f)*c( f k( f ) 



(36) 



~ (i) 

At lowest order, S L j is just the square of eikonal vertices, with elements given by 
color traces of Cj and c* L . 



4.2 Diagonalization of the Color Basis 

In lowest (linear) order in a s for the anomalous dimension matrices, the ordered 
exponentials of Eq. (]35[) may be reduced to sums of exponentials, by changing the 



color basis to one in which T s is diagonal. We proceed as follows. 

First, we treat the original basis of color structures as a set of kets, ||c^j}|. 

Second, we compute the anomalous dimension matrix T s at one loop in this basis 
and we then diagonalize it, finding its eigenvalues and eigenvectors. Let the diagonal 
basis of the eigenvectors be denoted as the set ||e^ K ) j, where from now on we shall 
use Greek indices for vectors and matrices expressed in this basis. The matrix defined 

as 

< C %| C <<> 3 ).(e%) x , (37) 
diagonalizes the anomalous dimension matrix according to the equation 



= A« V (38) 

Third, we reexpress the hard amplitudes in the eigenvector basis. In the previous 
notation, 



/><<>)• = (h<i>y (««>) 



t 



hf = {R<*) e hf. (39) 
We reexpress Eq. (^) in the new basis as 

T r (ff<')S<')} = K f >)*SS/4« = <sW, (40) 
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where we define 



o(f) 
D 7/9 



R 



(f) 

■yM 



] 9 {f) (r^ 

°MN \ n - 



In the diagonal basis of the eigenvectors, Eq. 



N/3 

I) then becomes 



(41) 



Trjtf( f ) 



' MjJ *V,*.tf)) S^{^,Ay,a s ^)) 



H$ ( Mil, Ay, a s (A S® (l, Ay, a s (l^j/N 1 



'/St 



// 



x exp 



m jj/ n dfi' 



(42) 



in which the leading iV-dependence of the product is organized in a sum of exponen- 
tials. 

In general, the eigenvalues in Eq. fl42f) depend on the flavors and directions of the 
colliding partons, both incoming and outgoing. In the following two sections, for the 
partonic processes f, we will give: 

• The basis of color structures ||c®/)|. Generally, we will work with an arbitrary 
number of colors, iV c . 

• In this basis, the anomalous dimension matrix T s . 



The set of eigenvalues, j Ajp j , and the corresponding set of eigenvectors, j | e^ K ) j , 
of Eg . Each eigenvector, in the notation of Eq. (|37|), corresponds to a column 
of^))" 1 . 



With these quantities, combined with explicit hard-scattering functions, it is possible 
to write down resummed jet cross sections f2~lf . 

The anomalous dimension matrices for the processes qq — * qq and qq — > qq have 
been known for some time |7], || [10| . In Sec. [5] we will review these matrices at one 
loop, and specify their eigenvectors. There we also give corresponding results for the 
process qq — > gg (related by time-reversal invariance to gg —* qq as well [17, [18]]) and 
for the process qg — > (related to — » gg). The quark-gluon processes involve a 
3x3 matrix, the quark-quark a 2 x 2, and all of these results are reasonably straight- 
forward. For gluon-gluon scattering on the other hand, the anomalous dimension 
matrix is 8 x 8 in 577(3), and requires a somewhat more extensive analysis, which we 
give in Sec. 6. 



5 Anomalous Dimensions for Quark-initiated Scat- 
tering 
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5.1 Formalism 



Our Wilson lines represent partonic processes, whose momenta and colors are labeled 
as in Eq. (pj), f A (I a, r A ) + Ib (Ib, r B ) -> fi (pi, n) + f 2 (p 2 , r 2 ). In terms of hadronic 
momenta, we take U = XiPi, i — A,B, with X{ the partonic momentum fraction. Re- 
ferring to the dijet process Eq. fll8| ) at partonic threshold, the dimensionless, lightlike 
velocity vectors t>f , vf = 0, which define the directions of Wilson lines at the vertices 
are given by 

It = MjX i — A,B 

$ = Mjjv? 2 = 1,2. (43) 

We also recall the definitions of partonic Mandelstam invariants, 



\2 



>2 



s = (l A + Ib) 
i = (I A - px 
u = {Ia~P2? ■ (44) 

The anomalous dimension matrices for all the processes below depend on logarithms 
of ratios of these invariants, and we introduce for them the following notation: 

T = ln(— 

s 

— ii 

U = ln( (45) 

s 

As mentioned above, although the full cross section is gauge independent, the 
functions into which it is factorized depend, in general, on the choice of n M , the axial 
gauge-fixing vector. From the factorized cross section, Eq. (p9|) , it is clear that gauge 
dependence in the product of and S^' must cancel the gauge dependence of the 
incoming and outgoing jets, ip and J^ l \ Because the jets are incoherent relative 

to the hard and soft functions, the gauge dependence of the anomalous dimension 

(f) 

matrices T s must be proportional to the identity matrix. To summarize this gauge 
dependence we introduce, for each parton in process f, the function 

with C/j = Cf (C a) for a quark (gluon), and with (see also Appendix 

_ (Vi ■ n) 2 



-- ImVj) In 2 H 

2 y 1 2 2 2 



(46) 



1 — 119 

\n\ z 

In these terms, we will write our anomalous dimension matrix as 



(47) 



(rf) KL = (tI))kl + | E %)(^))^l, (48) 

=A,B,1,2 
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ffl ({) 

and in the following, for specific processes, we will report only V s ;, to which T 3 
reduces under a proper choice of gauge. Eq. ( 48f) will allow us to recover the full 
gauge-dependence, whenever necessary. We emphasize here that the dependence of 
the anomalous dimension matrix elements on the external eikonal four-vectors and on 
the axial gauge-fixing vector comes entirely from the evaluation of the graphs shown 
in Fig. f|, no matter which partonic process we consider. The explicit computation 
of the graphs shown requires the eikonal Feynman rules, given in Appendix [A], and a 
few one-loop integrations, summarized in Appendix [B|. 

One last remark has to be made about the choice of physical channel s, t, or it, 
for the definition of the basis of color structures. All choices are allowed, and related 
to each other by simple crossing transformations. In the following, apart from the 
processes qq — > gg and gg — > qq, which are better described in terms of s-channel 
color structures, we will always use t-channel bases, which seems to be the natural 
choice when analyzing forward scattering M. In the following four subsections we 
present the results for the anomalous dimension matrices T s ; for partonic processes 
involving quarks. 



5.2 Soft anomalous dimension for qq —> qq 

We treat the process 

q {I A: r A ) + q (h, tb) — ► q (Pi,n) + q (P2, r 2 ) , (49) 
in the t-channel singlet-octet color basis 

c 2 — ~ 2N ^ TA,ri ^ rB < 7 " 1 2^ TA ' rB ^ ri ' T2 ' (50) 
The procedure we follow is the same as the one described in Refs. 0] and [17, 18] . We 



evaluate one-loop corrections for each color vertex, as shown in Fig. |], determining in 
every case the color decomposition of its ultraviolet divergences. The relevant scalar 
integrals are reviewed in Appendix [FJ below. The UV divergences found with the 
color vertex cf in Fig. |] in general generate counterterms for all the vertices c^p in 
the same set, where J may or may not equal I. This is standard operator mixing 
under renormalization. The elements of the anomalous dimension matrix are then 
found from Eq. (]34]) . 

In this manner, we find the matrix 



a q ( 2C F T 



71 




N, 



, Nc 

UT-2U) 



(51) 



The dependence on the logarithmic ratio T is diagonal in a t-channel color basis. 
Referring to the resummed iV-dependence in Eq. ([42]), we see that in the forward 
region of the partonic scattering (T — > — oo) color singlet exchange is exponentially 
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enhanced relative to color octet ||. The eigenvalues of this anomalous dimension 
matrix are 



{N 2 - 2)T + 2U - N C VK 



a s 1 

TT 2N C 

— ^W \( N c 2 - 2 ) T + 2U + N ^ 

7T Zl\ r 1 




(52) 
(53) 



Ai 
A 2 

where A is defined as 

A = N C 2 T 2 - 4TU + AU 2 . 
The corresponding (arbitrarily normalized) eigenvectors are 



ei 



(-2 



(54) 

From these eigenvectors we may reconstruct the matrix R' 1 of Eq. ([H]), which gives 
the transformation from the singlet-octet to the diagonal basis, for any scattering 
angle. The remaining anomalous dimensions are found in much the same fashion. 



5.3 Soft anomalous dimension for qq — > qq 

In this subsection we analyze the process 

q {I a, t a) + q (is, t b ) — > q (pi, n.) + q (p 2 , r 2 ) 
again in the t-channel singlet-octet color basis 

1 _ 1 

c i ~~ ~ 2N ^ TA ' ri ^ rB,r2 2^ rA ' T2 ^ rB ' ri 



We find the anomalous dimension matrix || 

-^(T + [/) + 2C F C/ 2U 



r.o/ 



7T 



■C7 



2CVT 



(55) 



(56) 



(57) 



Once more, the dependence on the logarithmic ratio T is diagonal in the t-channel 
color basis, and again the color singlet dominates the octet in the forward region. 
The eigenvalues of this anomalous dimension matrix are 



Ai 
A 2 



o.s 



Ti 2N C 

TT 2N r 



{N 2 —2)(T + U) — N C VA* 
(N 2 — 2)(T + U) + N c \fA! 



{hi 
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where A' is defined as 

A' = N c 2 (T-U) 2 + 4TU, 
and the corresponding eigenvectors are 



(59) 



ei 



-n c 3 (t-u)-2n c u-n c 2 va' 
(n c 2 -i)u 



e 2 



( -n c 3 {t-u)-2n c u+n c 2 ^a 
(n c 2 -i)u 

V 



(60) 



5.4 Soft anomalous dimension for qq —> gg and gg —> qq 

Next, consider the process 

q(lA,r A ) + q(lB,r B ) — > g(pi,n) + g(p2,r 2 ) , 
in the s-channel color basis 



(61) 



Cl 

c 2 
c 3 



jr 1 r 2 ccr r c\ 

■ F)r B r A 



d rir2C (T c F 

if rir2c (m rBrA , 



(62) 



where the Tp's are the generators of SU(N C ) in the fundamental representation, while 
jabc an( j ^abc are totally antisymmetric and symmetric SU(N C ) invariant tensors 
respectively. In this basis, we find the anomalous dimension matrix 



a. 



r« = - 



7T 






2 (U - T) 



U-T 
^(T + U) £a(U-T) 

%±(U-T) ^(T + U) 



(63) 



The same anomalous dimension describes also the time-reversed process p7| , [X8| 
9 (pi, r x ) + £ (p2, r 2 ) — ► q (l A , r A ) + q (l B , r B ) . 
The eigenvalues of IV are the solutions of the cubic equation 



(64) 



given by 



a 



^C A (T + U) X + [ — 



7T 



a. 



7T 



jV c 2 + 4 
4 



(U-Tf 



+ 



C A (T + [/)([/- T) 2 = 



(65) 



^2,3 



2±l( X W-Y + C A {T + U) 



Og_ 1 
7T 3 



_I( X V3 _ y) + c A (T + [/) ± ^(x 1 / 3 + y) 



(66) 
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where 



X 



Ca 
2 



3v/3 



(T + U)j +-C A (T + U)(U-T) 2 (N 2 c -8) 

6 (iV c 2 + 4) 3 _ {u _ t) 2{Ca {t + ^ 4 {N 2 _ 4) 



(U-T) 
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+ \ (U - T) 4 (T + U)) 2 ((iV 2 - 8) 2 - 12(iV c 2 - 2 



,1/2 



and 



Y 



im , tt\\ , 3 lTT m\2 / \-2 



-y(T + C/)l +-(U-TY(N' c +i) 



For each eigenvalue Aj the corresponding eigenvector is given as 



N C (U-T) 
2\ l -C A (T+U) 
1 



, 1 = 1,2,3. 



5.5 Soft anomalous dimension for — > and qg — > qg 

We now turn to the "Compton" process 

q (Ia, ta) + 9 (h, r B ) — ► q (pi, n) + g (p 2 , r 2 ) , 
in the t-channel color basis 



C2 = cf^ c (l|) n 
c 3 = ir^ c {T c F ) n 

The result of our calculation is 

/ (C F + C A )T 



TA' 



71 



V 



o 

2U 



C F T + ±fU 



2N C 



£a, 
2 

■u 



u 

2 u 

CpT + £±U 



which applies to the process 

q (pi, n) + g (pa, r 2 ) — ► g (/a, r A ) + g (l B , r B ) 



(67) 



(68) 



(69) 



(70) 



(71) 



(72) 



(73) 



as well. As in Eq. fl5l|) , the dependence on the logarithmic ratio T is diagonal in the 
t-channel color basis, and the ^-channel color singlet dominates in the forward region 
(T--00). 
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The eigenvalues have the same structure as in Eq. 
> a., 1 

•^2,3 



ifi (X' 1 / 3 -Y' + (3C F + C A )T + C A U 

7T 3 V 

" s 1 { -\{X'^ - y') + (3C F + C A )T + c A u ± ^(x' 1 / 3 + Y') 



7T 3 

with X' and Y' given by 



(74) 



X' 



3^ 



^ U 6 ( N ? + *) 3 _ U 2I Ca i t _"\\ (iV2_4) 
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2 



(cu (r - 1)) 2 ((x c 2 - s) 2 - i2(x 6 2 - 2; 



1/2 



and 



r = - 



(75) 
(76) 



The eigenvectors have a form similar to Eq. (|6 



u 



Xi-(C F +C A )T 

NgU 
2\i~2C F T-C A U 
1 



i = 1,2,3. 



(77) 



6 Soft anomalous dimensions for gg 
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In this section, we present the anomalous dimension matrix for the eikonal version of 
gluon-gluon scattering. Here the proper choice of a set of color vertices is somewhat 
less obvious. We would like to compute, as easily as possible, the color structure of 
one-loop corrections. This was not a problem for 2^2 partonic processes involving 
quarks, because the couplings of the fundamental representation are rather simple. 

In subsection |0| below, we introduce a basis of nine color tensors |2R |2D|, which, 
although overcomplete, is convenient for computation of the one-loop IV for the 
gluon-gluon case. We give a few details of how to color-decompose the one-loop 
graphs into this set, and present a first, "raw" version of the anomalous dimension 
matrix Ts>- In subsection |6.2| , we modify the original basis, to one in which IV is in a 
block-diagonal form in which one block is diagonalized. We also present eigenvalues 
and eigenvectors in this basis. Finally, in subsection |6.3| , working with N c = 3, we 
recall t-channel color projectors describing the decomposition of the product of two 
color octets into irreducible representations |30fl . We use these projectors to reduce 
our set of color structures to a complete, eight-dimensional basis of color tensors. We 
exhibit the resulting one-loop anomalous dimension matrix, along with its eigenvalues 
and eigenvectors. As in all cases above, logarithms of t/s appear only in the diagonal 
of IV m this basis. 
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6.1 Initial basis and Ts> at one loop 



To begin with, we consider the set of color tensors, in terms of which an arbitrary 
four-gluon diagram, with external color indices r^, rg, r\ and r 2 can be expanded 



29] 



c x = Ti(T r F A T r F B TpTP) 

c 2 = Ti(T r F A T F B T F 'Tp) 

c 3 = Tr (T F A T F T r F T r F B ) 

c A = Tr (T F A T F T r F T T F ) 

c 5 = Tr (T F A TpT r F 'T r F B ) 

c 6 = Tr (T F A T F 2 T r F B T F ) 

C9 A r2^r B ri 1 



(7f 



where the Tp's are the generators of 5*f/ (iV c ) in the fundamental representation. This 
set is illustrated graphically in Fig. ||], where all the lines and vertices have only color 
content. The correspondence between the picture and Eq. ( |78[) is straightforward, in 
terms of the color dependence of the quark-gluon vertex in the QCD Lagrangian. 

As in the case of scattering represented by Wilson lines in the fundamental rep- 
resentation, there is no mixing of operators with different numbers of Wilson lines. 
Therefore we only need to consider color tensors with four indices in the adjoint rep- 
resentation. Eq. (|78|) clearly includes the maximum number of inequivalent traces we 
can build out of four generators T F , in addition to the three possible singlet com- 
binations. The set (78) mixes only with itself under renormalization. To see this, 
consider one- loop corrections to the vertices of Eq. ([78]). The color decomposition 



of these one-loop corrections is obtained by the graphical identities shown in Fig. 
whose analytic form is 



r 

E 

k 



T2T3 



-2i [Tr (TpTpTp) - Tr (TpTpTp)] 



T4T2 



2N r 



J rsri w r4T2 i 



(79) 



where f rir2T3 are the structure constants of the SU(N C ) algebra. The first identity 
allows us to replace three-gluon couplings by color traces in the fundamental rep- 
resentation, and the second to "expand" the color content of internal gluons and 
Wilson lines into pairs of color lines in the fundamental representation. External 
gluon lines may be treated in the same way, by use of the identity 5 a b = 2Tr T F T F . 
It is straightforward to extend this procedure iteratively to an arbitrary order, using 
the corresponding identities for the four-point function. Essentially the same color 



decomposition can be used to show that the basis ( 78f) is complete in the expansion 
of the hard-scattering function for gluon-gluon scattering, as in Eq. (E32|) above. 
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As in Sec. 5, we follow the procedure of Refs. M and |T7], p 



and compute the 
8|) of the one-loop 
3) is 



anomalous dimensions from the color decomposition into the set 
ultraviolet divergences of Fig. [|, using Eq. fl3"3|) . The choice of basis in Eq. 
particularly well-adapted to this procedure, when we employ the identities of Eq. 
([79D. After straightforward, although rather tedious, calculations, we find the 9x9 
anomalous dimension matrix, 



7T 



T 






T-U 

N c 

_U_ 

N c 

V o 





u 







_T_ 

N c 
U-T 





T 




T-U 

N c 

U_ 

N c 







(T + U) 



T_ 

N c 



2- 

Nc 








u 





— 2L 

N c 
U-T 
N c 











N c 



u_ 

N c 



u_ 

' Nc 



U_ 

N c 

U_ 

Nc 







(T + U) I 



2T 





T-U 

Nc 
U-T 

N c 

T-U 

N c 



U-T 

N c 








_T_ 

N c 



T_ 

N c 
_T_ 
N c 
T_ 

N c 




2U j 



Interestingly, this anomalous dimension matrix is diagonal in the large N c limit. 



(80) 



6.2 Block-diagonal form of the anomalous dimension matrix: 
eigenvalues and eigenvectors 

Although it has many zeros, the anomalous dimension matrix of Eq. (^) is not 
transparently easy to diagonalize for arbitrary N c . As a first step in its simplification, 
we will make a change of color basis that transforms (^) into a block-diagonal form. 
In particular, looking at the first three rows of the matrix, we notice that, if we are 
able to eliminate the non-zero off-diagonal entries, we immediately determine three 
of the eigenvalues. 

The way to accomplish our purpose is to employ symmetric and antisymmetric 
linear transformations from the set (|78| ) to a new set {c' 7 }, 





= Tr (Tp A T r F B Tp Tp 1 ) 


-TiiT^Tp'TpTp 3 ) 


= Ci — c 3 , 


/ 

C 2 


= Tr(T F A T r F B Tp Tp ) 


-Ti(T F A Tp 2 Tp l Tp B ) 


= c 2 - c 5 , 


c's 


= Tr (Tp A TpT F B Tp 2 ) 


- Tr (Tp A T r p 2 Tp B T r p 1 ) 


= C4 — c 6 , 


i 

c i 


= Tr(r; A T^TpTp 1 ) 


+ Tr(T r p A Tp 1 TPT r p B ) 


= Ci + c 3 , 


i 

C 5 


= Tr (T F A T F B Tp T T p ) 


+ Tr (T F A Tp Tp T F B ) 


= C2 + C5 , 




= Tr (T F A Tp 1 T F B T F 2 ) 


+ Tr (Tp 4 T F 2 T F B Tp ) 


= C4 + Cq , 


/ 

C 7 


_ 1 r r _ 

A^rA,ri^ r B r2 c 7 j 






c's 
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Cg. 



The expression of IV in this new set of color structures has the form 

rv = 

where the block r 3x3 is diagonal, 



1^3x3 3x6 

06x3 r 6x6 



(82) 



T - tts 

1 3x3 — 

7T 









NJU 








V N C (T + U) 



while the block r 6x6 is given by 



6x6 





/ N C T 








-[/ 


T-U 


o \ 







N C U 








U-T 


— T 










N C (U + T) 


[/ 





T 


7T 


2 (T - U) 





2T 


2N C T 










-2U 


-2T 
















V o 


-2 (T - U) 


2C7 








2N C U J 



?4) 



In both the previous matrices we have used explicitly = iV c (compare with (|80|)). 

The first three eigenvalues, Ai, A 2 and A 3 , can be read directly from Eq. (p3|). To 
determine the others, we just need the diagonalization of a 6 x 6 matrix. We have 
computed its eigenvalues and obtained 



Ai 
A 5 
A 6 
A 7 

A8,9 



-N C T = A x , 

7T 

-N C U = A 2 , 

7T 

— N c (T + [/ 

7T 
Os 

7T 
«s 

7T 



A 3 , 
2 



X l/3 _ y + Nc (T + [/) 



( X V3 _ y) + at c (T + C/) ± -iV3 (X 1 / 3 + F 



$5) 



where, as before, to reduce clutter in the notation, we introduce auxiliary quantities 
X and Y, defined for this particular 



x 



-N e (Nl 
4 



9 



2 (T 3 + f/ : 



3{T 2 U + TU 2 



9 



12 (X 2 - l) [T 6 + f/ 6 - 3 (T 5 U + TU l 



-3 ( N? + 6iV,? + 33JV? + 24 ) ( T*U' Z + T Z U 



1/2 



.4rr2 



n2rr4 



+6T 3 U 3 (n® - 4X C 4 + 53N 2 + 14) } 



56) 



24 



and 



Y = --(N* + 3 



T + U — TU 



x- l '\ 



(87) 



The eigenvectors corresponding to the first three eigenvalues have the form 



e {3) 

(6) 



, i = l,2,3, 



e< 3) = 



(3) 

where the superscripts refer to the dimension. Thus, the e\ s are three-dimensional 
vectors, defined as 

( 6a \ 

S i2 , % = 1, 2, 3 , (89) 

V <*« / 

while O^ 6 -* is the six- dimensional null vector. 

The eigenvectors corresponding to the other eigenvalues have the general form 



4,. ..,9. 



(90) 



In particular ef\ ef 1 and ef^ are given by 



e 4 



/ (jV e 2 -2)T \ 
2U 

1 

T-U 



N, 



u 

U-T 



\ 



u 
-N c 





p (6) 
e 5 



V 



{N*-2)U 
2N C T 

T-U 

TN C 



1 

U-T 



p (6) 
e 6 



T 

' U 



) 



V 



[Nl - 2) 


-N r 



U-T 
2U 



(91) 



while ef\ eg 6 " 1 and ef^ are given in terms of the coefficients of a s /7r in the correspond- 
ing eigenvalues, 

(92) 



7T 



K — \ , 



by 



7,8,9. 



(93) 



/ *i(AJ) \ 
,(6) _ x 3 (\$ 

1 

V z 6 (a;) y 

Here, once more to shorten our rather complicated formulas, we have parametrized 
the eigenvectors in terms of the quantities 

(\/\ 1 f y | r 4iV c (T 2 - U_) - 2AK2T - tQ ] 

UJ " 2C/\ l ^ [N c (T + U)- X'i\ (X'i - 2N C U) J ' 



25 



x 2 (K) 

x 6 (A-) 



2T 
1 



4UT 2 



i\> + 2Ux 1 (\>)} 1 

{[K + 2Ux 1 (\>)} [(AJ - 2N C U) (K - N C U) — 2T (T — 17)] 



1 
T 



2T (T-U) (K-2N C U)}, 
^[(N c T-X' l )x 1 (X' i )+T-U], 
(K - N C U) [A - + 2Uxi(X' i )] 



2T 



+ U-T 



Working with N c = 3, the eigenvalues simplify to 
Ai 



A 2 

A 3 
A 7 

A8,9 



A 4 = 3^T, 

71 

A 5 = 

7T 

A 6 = 3^(T + f7), 

7T 

2— (T+U), 

7T 

"* r 2T + 2U ± AVT 2 + U 2 -TU 



7T 



For N c = 3, the eigenvectors , and eg '' become (compare Eq. (pi]) ) 



(6) 



J6) 
e 4 



/ _7T \ 
' 6£7 » 

_ 1 

3 

U-T 
3U 

T-U 
U 

1 



1 T 



,(«) 



3 T-U 

7U 
6(T-U) 

_ 1 

3 



T-U 



\ 



,(6) 





/ 


V 


1 / 


while e\ 6 , eg 6 ' 1 and eg 6 ^ reduce to 












/ ai (AQ \ 




-1 




a 2 (A' i ) 


p (6) - 


-1 






e 7 — 


1 
1 

V i ) 




a 4 (A' J ) 
1 

V o 6 (AD J 



T 
3U 
_ 1 
3 

7{T-U) 
6U 
_T 
U 



1 



8,9. 



Here for brevity we have introduced the quantities 

-1 
6UK 

ai(Aj,T<-I7) 



or (AD 
a 2 (A-) 



20(T 3 - U 3 ) - 52TU(T - U) + (10T 2 + 5f7 2 - 9TU)X[ 
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a 3 (X) = — \— {-36TU(T 2 + U 2 ) + 32T 2 U 2 + 20(T 4 + U 4 ) 
bl U K 

+ [10(T 3 + U 3 ) + 5TU{T + £7)]AJ}, 

a 4 (AO = — ^— [lOf/ 4 + 20T 4 - 54T 3 [7 + 58T 2 [/ 2 - 34Tf/ 3 

+ (10T 2 + 7U 2 — 7TU)TXj\, 
a 6 (Aj) = a 4 (A<,T<-l7), (9? 



where the A'/s are obtained from the Aj's of fl95|) according to Eq. ( [92]) . The quantity 
if is defined as 

K = 5(T 2 + U 2 ) - 6UT. (99) 

Using the formalism of subsection ^]2|, we can write the matrix , having the 

eigenvectors of the anomalous dimension matrix in its columns. We have verified that 

is a diagonal matrix with the eigenvalues of Eq. fl95|) . 
6.3 Color projections for gg —> gg 

So far we have been working with the sets of color structures in Eqs. Q78| ) and (Wfh. 
These sets are both overcomplete, since the c/'s are not linearly independent ||28|| , 
as we will show below. In this section we would like to rewrite the set in Eq. flgip 
in terms of SU(3) tensors (we keep N c = 3), and to elucidate the group structure 
of t-channel color exchange for the product of Wilson lines that generates the same 
noncollinear soft radiation as gg — > gg. 

To relate the basis of traces to color exchange, we first recall the product formula 
for generators of SU(N C ) in the fundamental representation, 

t f t f = U l3 \ + 1 (d ijk + if ijk ) T k F . (ioo) 

Eq. ( |100| ) enables us to rewrite the set of color structures in Eq. fl8T|) in terms of the 
tensors / and d, 

c i = Tr (T F A T F B T F 2 T F ) — Tr (T F A T F T F T F B ) = — [fr A r B id ri r 2 i ~ d rArB if rir2 i] , 
c 2 = Tr (T F A T F B T F T F ) — Tr (T F A T F T F T F B ) = — [f rA r B id ri r2i + d rArB if rir2 i} , 
c'3 = Tr (T r F A T r F T F B T r F ) — Tr (T F A T F 2 T F B T F ) = - [f rAri id rBr2 i + d rAri if rBr2 i] , 

C4 = Tr (T F A T F B T F 2 T F ) + Tr (T F A T F T F 2 T F B ) = ~Sr A ri^r B r2 T [dr A rildr B r 2 l 

o 4 

fr A r\lfr B r 2 l\ j 

C 5 = Tr (T F A T F B T F T F 2 ) + Tr (T F A T F T F T F B ) = ~^^r A r B ^r 1 r2 ~\~ T [dr A r B ld ri r 2 l 

o ' 4 

fr A r B lfrir2l\ j 

c% — Tr (T F A T F 1 T F B Tp) + Tr (T F A T F 2 T F B T F 1 ) = —S rAn 8 rBr2 + - [d rAri id rBr2 i 

o ' 4 
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fr A rilfr B r2l] j 

4' 

, I 

C g = -^r A r 2 ^r B ri • (101) 

We now go on to discuss color in the t-channel. 

The color content of a set of two gluons in SU (3) is described by the direct product 
8 ® 8, which can be decomposed into irreducible representations in the standard way, 
as 

808 = 1 + 85 + 8^ + 10 + 10 + 27. (102) 



In Ref. [3(| we can find the representation in terms of SU(3) tensors of the projectors 
performing the above decomposition. In our notation, t-channel projectors are given 
by 

Pi(r A ,r B ;r 1 ,r 2 ) = -S rAn Sr B r 2 , 

3 

P8 s ( r A,r B ;r 1 ,r 2 ) = -d T ATxC d r Br , zC , 



P8 A ( r A, fB] f\i ^2) — ^fr A ricfr B r 2 c, 



^"l0+10 vAi r Bi r l 1 ^2) 2 \^TATB^riT2 ^r A r 2 ^r B ri) ^frATxcfr B r'zci 



PwiT Ai r B\ r \-, r 2) — ^;(^r A r B ^r 1 r2 + ^r A r2^r B n) — Q^r A ri^r B r2 



1 

2 l 
3 

~zd r A ncdr B r2C • (103) 

5 



From a comparison of Eq. ( |101| ) and Eq. ( |103| ), we see that the first three elements 
of the basis ( 101|) , depending only on the "mixed" combinations fd, df, cannot be 
expressed in terms of these projectors. On the other hand, it is easy to find relations 
for the remaining basis tensors, 

c' 5 = --Pa - -P 8s + -P 27 , 
C ' 6 = 3 Pl + 12 P8s " 4 P8a ' 



c' 7 = 2P 1 
1 

4 



c's - - fPi + P8 A + Ps s + PioeTo + P27 J • 



c'. 



\ (A - flu + Ps s - PioeTo + P27) ■ (104) 
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Evidently, the six tensors d T , i = 4 ... 9 may be replaced by only five t-channel 
projectors. This is because, as pointed out above, the original basis of Eq. (|78|) is 
overcomplete, and indeed, the basis elements of Eq. fl31~| ) satisfy the relation 



1,1,1 

C4 + C 5 + Cq = C 7 + C S + Cc 



(105) 



In terms of the invariant tensor d^, related to the c'j in Eq. ( |101|) , this is equivalent 
to the identity |j§ 



jirn^irnk 



- (5ij5ke + 5jk5a + 5ik5ji) 



(106) 



It may be worth noting that the Jacobi identity, as well as a related identity involving 
products of /'s and cfs in Ref. is simply an expression of the cyclicity of traces 
in the original basis, Eq. (f78p, and hence does not result in a further reduction of its 
dimension. 

Our new basis is thus given by the eight color structures 



jc'i, c' 2 , c' 3 , Pi, P$ s , P$ A , -P 10 eTo) -^27} 
In this basis, the soft anomalous dimension matrix becomes 



r.«/ 



1^3x3 03x5 
05x3 T5x5 



(107) 



(108) 



with r 3X 3 given by Eq. (p3|) with N c = 3, and with r 5X 5 given by 



5x5 



7T 



V 



6T 





-QU 











3T + ^ 


— EL 


-3U 





3U 
4 


3U 
2 







9f7 
4 





6U 
5 





3U 


9U 
5 








2U 
3 


AU 
3 


-2T + 



(109) 



From Eqs. (|83|), (|108|) and ( |109| ) , we see that the anomalous dimension matrix is now 
diagonal in its T dependence. 

We have computed the eigenvalues of the matrix above, to check that they are 
the same as in Eq. (|95|) (we will have one fewer, due to the dimensional reduction of 
our matrix), finding, indeed, 



A 4 
A 5 
A 6 
A 7 
A s 



A 1 = 3^T. 

7T 

A 2 = 3^U 



A 3 = 

a s 
2— 

7T 

2— 

7T 



7T 

3—(T + U) 

71 



T + U - 2VT 2 -TU + U 2 



T + U + 2VT 2 -TU + U 2 



(110) 
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The eigenvectors have the general form 



(5) 



1,2,3, 



where the e- 3 '* are given explicitly in Eq. fl39|); ef\ ef* and ef } are given by 



(3) 

e (5) 

(5) J5) 



4...; 



(Ill) 



,(5) 



p (5) 
e 4 



/ -15 \ 

u 2 U 
_ 15 T 

2 1/ 

3 
1 



V 



/ 



,(5) 



3 _ 3 T 

4 2U 



V 



,(5) 



15 \ 



3 _i_ 15 T 
' 2 *~ 2 U 



15 
2 



15 T 
2 U 



\ 



-3 
1 



(112) 



ef^ and eg° j have a rather complicated form, which can be simplified by the reparam- 
eterization, 



,(5) 



,(5) 



( HK) \ 

b2(K) 

HK) 
l 



7,; 



(113) 



where A'j was defined in Eq. 

3 , 



UK) 
UK) 

HK) 



U 2 K' 



and where the b^s are given by 
80T 4 + 103C/ 4 - 280UT 3 - 300TU 3 + 404T 2 £/ 2 



+ (40T 3 - 16f/ 3 - mT 2 U + 52TC/ 2 )A' i ] , 



2K 



7 [20T 2 - 50C/T + UU 2 + (10T - 5C/)A' i ] 



-[40T 3 - 64t/ 3 - 120T 2 t/ + 130TC/ 2 



2C/K' 

+(20T 2 + UU 2 - 20TC/)A-] , 
— (2T + 5C/-2A',), 



-ft'' being the auxiliary function 

K' = 20T 2 - 20UT + 21U 2 . 



(114) 



lis) 



As T — > — oo, the eigenvectors e,i approach simple linear combinations of the color 
tensors of Eq. ( |107| ). In particular, becomes proportional to the color singlet vector 
(0, 0, 0, 1, 0, 0, 0, 0) in Eq. ( |107[ ), with an eigenvalue that decreases as 6T. This is the 
most negative eigenvalue, so that, once again, color singlet exchange dominates in 
the product of hard and soft functions, Eq. (f|^), for gluon-gluon scattering. Of the 
remaining eigenvalues, the color structures of the first three are simple products of 
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the / and d tensors for all angles (see Eq. ( |101| )). Next, from Eq. ( 112 ), we see 
that e4 and e§ become degenerate and approach linear combinations of symmetric 
and antisymmetric octet projectors, while corresponds to 10 © 10 exchange in 
the forward limit. Finally, eg approaches pure 27 exchange in the forward limit, 
with the largest, and hence most suppressed, eigenvalue, — 2T. In fact, suppression 
increases with the dimension of the color representation exchanged. This suggests 
that radiation by accelerated charges in a nonabelian gauge theory increases with the 
dimension of the color exchange, not only from singlet to octet, but also to higher 
representations. 



7 Conclusions 

The anomalous dimension matrix for soft radiation in gluon-gluon scattering com- 
pletes the set of one-loop calculations that we set out to analyze. An interesting 
regularity of our results, relevant to near-forward scattering, is that for each flavor 
combination, logarithms of t/s appear only in diagonal matrix elements in those bases 
that describe definite color exchange in the t-channel. This had been shown previ- 
ously for the cases of quark-quark and quark-antiquark scattering at one loop || [| 
and two loops flOfl . One application of the anomalous dimension matrices and their 



eigenvalues is for the resummation of threshold singularities to next-to-leading loga- 
rithm in jet cross sections |21] at fixed rapidities. This was the original motivation 



for our study. We are hopeful, however, that the formalism of resumming coherent 
noncollinear soft gluon radiation in terms of color mixing will have a variety of other 
uses as well. 
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A Feynman rules for eikonal diagrams in axial gauge 

In this Appendix we summarize the Feynman rules we have used in the evaluation of 
the eikonal diagrams of Fig. We refer to Fig. [7|. 

The propagator for a quark, antiquark or gluon eikonal line (Fig. 0) is given by 



Sv ■ q + ie 

where 5 = +1(— 1) for the momentum q flowing in the same (opposite) direction as 
the dimensionless vector v. The interaction vertex for a quark or antiquark eikonal 
line (Fig. 0a, b) is 

-ig s T c Fba Av», (117) 
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with A = +1(— 1) for a quark (antiquark). 

For the gluon eikonal vertex (Fig. |7|c), we may take 



-g s f abc Av„ (118) 

where we agree to read the color indexes a, b, c anticlockwise, and where A = +1(— 1) 
for the gluon located below (above) the eikonal line. 

Finally in our calculations we employ the general axial gauge gluon propagator, 

D^{k) = —^-N^{k) J N*»(k) = <r- , +^7^4^ (H9) 

At + ie n ■ k (n • ky 

with the axial gauge-fixing vector. 



B Evaluation of one-loop eikonal vertex correc- 
tions 

We refer to Fig. £|, where we have represented the UV divergent 0(a s ) contributions 

ff] (f) 
to the eikonal vertex Wj . The counterterms for W} are ultraviolet divergent co- 
efficients, depending on the external momenta and on the axial gauge fixing vector, 
times our basis color tensors, c) . 

At Born level the eikonal vertices describe only color flow, and the relation 

WfLn = c? (120) 

holds. 

The one-loop corrections can be written in matrix notation as 

w$£ = <M*z 8 = w$£ Zgt (i2i) 

where the superscript T indicates transpose. The dimension of the vectors (and 
therefore the rank of the matrix of renormalization constants) depends on the specific 
partonic process, as discussed in detail above. 

We denote the kinematic part of the one- loop vertex correction to Wj , with the 
virtual gluon linking lines V{ and Vj, as Uij(5iVi, SjVj, Aj, Aj), with Si's and Aj's defined 
in Appendix |A|. Its expression is 



^(SiV^SjVj.A^Aj) = gl 



2 , d n q —i [ Aj Aj Vi-Vj 



(2n) n q 2 + ie [ (5iVi-q + ie)(5jVj-q + ie) 

Ajfj-n P AjVj-n P P 

' + n 



(8iVi-q + ie) (n-q) {5jV r q + ie) (n-q) (n-q) 2 J ' 

(122) 

where P stands for principal value, 

, P x = - (- - T + (-l) /3 7 -,V (123) 

(q-ny 2\(q-n + iey K ' {-q-n + ief V ' 
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We rewrite ( 122 ) as 



1 1 

h{8iVi, SjVj) - -I 2 (SiVi, n) - -I 2 (SiVi, -n) 



11 

-^h^jVjin) - -hiSjVj, -n) + h(n 2 



where S^ is an overall sign 

= Ai Aj 5i 5j. 
The ultraviolet poles of the integrals are [0, § 



UV pole 



jUV pole 
L 2 



jUV pole 
J 3 



tUV pole 
J 4 



m s r 2 l 

e 
1 



7T If 



7 + In 5i5j 



Vi ■ Vj 



ln(47r) 



0:.s 



2vr U 2 e 
2 1 



2tt U 2 e 
7r e 



[7 + ln(^) - ln(47r)] 
[7 + ln(^) - ln(47r)] 



where all the gauge dependence is through the variable 

(vi-n) 2 



\n\ 



The double poles cancel in ( |124| ), giving 



a-.. 



Uij(SiVi,6jVj,Ai,Aj) = -Sij — 

7T6 



5i 5j v^Vj 1 
l n ( — ^_ — i) _ _ ln(^^) + 1 



(124) 
(125) 



(126) 



(127) 



'128) 



In order to obtain contributions to the different entries of the matrix of renormaliza- 
tion constants, the above expression has still to be multiplied by the color decompo- 
sition of its corresponding diagram into the basis color structures [0, [H], [T8| . 
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A 
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Figure 2: Identity representing factorization in leading regions of corrections to the 
eikonal vertex cj. S represents lines of zero momentum, while A, B, 1 and 2 label jet 
sub diagrams. 
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Figure 4: Eikonal vertex corrections for partonic processes contributing to the soft 
anomalous dimension matrices. 
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Figure 5: Graphical representation of the initial color basis for the gg — > gg process. 
All lines have color content only. 
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Figure 6: Graphical representation for the color decompositions in Eq. (|75|). 
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Figure 7: Illustration of the eikonal Feynman rules for a quark (a), antiquark (b) and 
gluon (c) eikonal line. 
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